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of cosmological singularities can be realized. Amongst others, it is possible to describe 
inflation with the appearance of a Type IV singularity, and phantom late-time acceleration 
which ends in a Big Rip. Finally, for completeness, we also present the Type IV realization 
in the context of suitably reconstructed F(R) gravity. 
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1 Introduction 

The observational data concerning the effective equation-of-state (EoS) parameter of dark 
energy w e g, slightly indicate that in the recent cosmological past the EoS parameter might 
have crossed the phantom divide and lies between quintessence and phantom regimes [1-4] , 
and therefore the cosmological evolution of the present Universe might be described by an 
effective phantom phase. There might also be the possibility that although the current 
Universe may be of quintessential type, the phantom era may occur in the near future. For 
an incomplete list of articles discussing the possibility that our Universe is described by 
an effective late-time phantom era see [1-16] and references therein. On the other hand, 
the inflationary era is also an accelerating phase of the Universe, which is not up to date 
excluded to be of quintessence or phantom type. One of the biggest problems in this type 
of phenomenology, is to fully understand the relation and the connection of early-time 
inflation with late-time acceleration, and in particular how to consistently describe these 
accelerating eras using the same theoretical framework. The origin of the problem itself 
might be the incomplete understanding we have, up to date, of dark energy and how the 
transition from deceleration to the dark energy acceleration epoch is eventually achieved. 
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As was demonstrated in [1], a promising theoretical unified description of both early- 
time phantom inflation and late-time phantom acceleration is achieved by using canonical 
and non-canonical scalar fields. Using this generalized multiple-scalar-field scalar-tensor 
theory, with generalized scalar-field-dependent kinetic functions in front of the kinetic terms 
and with a generalized scalar potential, a consistent theoretical framework can be provided 
for these theories [17]. This can be done either by using two scalar fields, or by using even a 
single scalar field. However, in the latter case there is a theoretical obstacle connected with 
the transition time at which the EoS evolves from non-phantom to phantom and vice-versa. 
In particular, the single scalar field theoretical description of cosmological dynamics always 
comes with the cost of possible instabilities at the cosmic time that the phantom divide 
is crossed, since at exactly this point the dynamical system that describes the evolution 
of the Universe can be unstable towards linear perturbations [18], rendering the resulting 
scalar-tensor theory unstable. Hence, the use of two scalar fields is necessary in order to 
resolve these issues in a concrete way [1, 19]. 

Definitely, the single-field scalar-tensor theory can in some cases prove valuable, since 
it offers the opportunity to describe exotic cosmological scenarios, which were not possible 
to be realized in the context of standard Einstein-Hilbert general relativity. In a recent 
study [20] we addressed the issue of having an inflationary era in the presence of finite-time 
singularities, with special emphasis in the Type IV singularity. Particularly, as we explicitly 
have shown, it is possible to have viable phenomenology even in the presence of Type IV 
singularities, and in some cases compatible with observational data coming from Planck 
[21] and BICEP2 [22] collaborations. These finite-time singularities can be of various types 
and were classified in [23], with regards to their relevance in the cosmological evolution. 
These types vary, with the main criterion that classifies them being the fact that certain 
observable quantities may become infinite at the time these occur. 

In principle, a singularity in a cosmological context is an unwanted feature of cosmolog¬ 
ical evolution. Space-time singularities were studied long ago by Hawking and Penrose [24], 
and they are defined as isolated points in space-time at which the space-time is geodesically 
incomplete. This feature practically means that there exist null and time-like geodesics, 
which cannot be continuously extended to arbitrary values of their parameters. The most 
characteristic example of this type of singularity is the initial singularity of our Universe, 
which naturally appears in inflationary theories [25-28]. It is worthy to mention however 
that there exist theories which are free of this initial singularity, such as the bouncing 
Universe scenarios [29-50]. Apart from these severe singularities however, there exist sin¬ 
gularities at which some observable quantities may become infinite, without this implying 
geodesic incompleteness. These are called sudden singularities and were extensively stud¬ 
ied by Barrow [51] (see also [52-74] for some later studies). In these singularities, the 
scale factor remains finite at the singular point, and therefore these are not singularities of 
crushing type, such as the Big Rip [10-12, 23, 75-92], 

In our recent study [20], the focus was on the existence of a non-crushing type of 
singularity, namely the Type IV. Our theoretical description of the cosmological evolution 
involved a scalar-tensor theory with a single scalar field. However, when someone uses 
a single scalar field, the scalar-tensor description might be unstable at the cosmic time 
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at which the phantom divide is crossed. This was not the case in our previous article, 
but this is a clear possibility when someone deals with single-field scalar-tensor theories. 
Therefore, the purpose of the present article is two-fold. Firstly, we aim to highlight the 
problem of instability in the single-field scalar-tensor cosmologies, in the presence of finite 
time singularities. Secondly, we are interested to formally address the crossing from the 
non-phantom to the phantom era consistently. The latter issue is definitely connected to 
the first one. We shall exemplify the problem by using a single-field scalar-tensor theory 
to produce a variant of a well known potential from inflationary cosmology, namely the 
hilltop potential [93-100], in which we consistently incorporate a Type IV singularity. As 
we show, the solution is unstable at the time the EoS crosses the phantom divide, and 
consequently the two-scalar-field cosmological description is rendered compelling. Having 
done this, we shall study illustrative examples of two-field scalar-tensor cosmologies, in 
which finite-time singularities exist. More importantly, we shall explicitly demonstrate that 
it is possible to unify phantom (or non-phantom) inflation, with phantom (or non-phantom) 
late-time acceleration. Finally, for completeness, we shall investigate which F(R ) gravity 
can generate a cosmological evolution corresponding to the hilltop scalar-tensor model, 
near the Type IV singularity. Interestingly enough, if the Type IV singularity is assumed 
to occur at late times, the resulting F(R ) gravity is of the form R+ a well known model 
that can consistently describe late-time acceleration [101]. 

This paper is organized as follows: In section 2 we show that the single-field scalar- 
tensor theory may lead to instabilities at the phantom-divide crossing, which are removed 
with the use of a second field. In section 3 we proceed to the investigation of various specific 
cosmological evolutions in the framework of two-field scalar-tensor theory, which prove the 
capabilities of the construction. In particular, in subsection 3.1 we reconstruct scenarios 
which exhibit Type II or Type IV singularities, in subsection 3.2 we analyze the slow-roll 
inflationary realization in such theories, in 3.3 we reconstruct a universe with a singular 
inflationary phase that results in a dark energy epoch which ends in a Big Rip, in 3.4 we 
present the two-field Barrow’s model, while in subsection 3.5 we extend our investigation in 
the case of multiple scalar fields. In section 4, for completeness, we also present the Type 
IV realization in the context of suitably reconstructed F(R) gravity. Finally, in section 5 
we summarize our results. 

2 Scalar-tensor gravity: general analysis 

In this section we shall examine in detail how the finite-time singularities can be consistently 
incorporated in generic scalar-tensor theories containing two scalar fields. In addition, we 
will demonstrate why the need for two scalars is in some cases compelling. Using well 
known reconstruction methods [1, 102], for a given cosmological evolution we shall provide 
the details of the generic scalar-tensor models that can consistently incorporate finite-time 
singularities, and then in the next section we will explicitly demonstrate our results using 
some illustrative examples. We mention to the reader that apart from Refs. [1, 102], many 
reconstruction techniques had been developed earlier, and most of these studies are very 
well reviewed in Ref. [103]. 
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In this work we consider a spatially flat Friedmann-Robertson-Walker (FRW) metric 
of the form 

ds 2 = —dt 2 + a(t) 2 ^2 (d% 1 ) 2 , (2-1) 

i= 1,2,3 

with a(t) the scale factor. In this case, the effective energy density p e s and the effective 
pressure p e g are defined as 

PeS = ^2 H 2 , PeS = ~~^2 ( 2 ^ + 3 H 2 ^j , (2.2) 

where H = a/a denotes the Hubble parameter and dots indicate differentiation with respect 
to the cosmic time t. 

Let us first recall the finite-time cosmological singularities classification. For details on 
this subject the reader is referred to [23, 104-106]. According to [23, 106], the finite-time 
future singularities are classified in the following way: 

• Type I (“Big Rip”) : As t —>• t s , the scale factor a, the effective energy density p c Q, 
and the effective pressure p e g diverge, namely a —>• oo, p e g —>■ oo, and |p e ff| —>• oo. 
For a presentation of this type of singularity, see Ref. [75] and also Ref. [23]. 

• Type II (“sudden”): As t —>• t s , both the scale factor and the effective energy density 
are finite, that is a —>• a s , p e s —>• p s , but the effective pressure diverges, namely 
\p eS \ -t oo [51-73]. 

• Type III: As t —>• t s , the scale factor is finite, a —> a s but the effective energy density 
and the effective pressure diverge, p e g —>• oo, |p e ff| —> oo [23]. 

• Type IV: As t —>• t s , the scale factor, the effective energy density, and the effective 
pressure are finite, namely a —» a s , p e s —>• p s , |p e fr| —> Ps, but the Hubble’s rate higher 
derivatives diverge, that is H = a/a [23]. 

Finally, note that in the aforementioned list of singularities, once should add the stan¬ 
dard case of recollapse due to a positive spatial curvature ending in the “Big Crunch” 
[107]. 

In a previous work [20] we studied how singular inflation can be properly accommo¬ 
dated to a power-law model of inflation, studied by Barrow in [74], Using the scalar field 
reconstruction method [1, 102], we found which scalar-tensor gravities can lead to a suc¬ 
cessful description of such inflationary cosmology. However, the conventions we used in the 
previous article lead to a stable scalar-tensor solution, which had an effective equation of 
state characteristic of quintessential acceleration, as long as the power-law inflation model 
is considered. Nevertheless, the stability of the scalar-tensor solutions is not ensured in 
general, and instabilities might occur during the reconstruction process. Therefore, in this 
section we explicitly demonstrate how instabilities can indeed appear, by using two phe¬ 
nomenologically appealing inflationary models. As we evince, the simplest way to avoid 
the instabilities is to use two (or more) scalar fields, with one of them being ghost and the 
other scalar field being canonical. 
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In a first subsection we aim to show how the Type IV singularities can occur in a 
slightly deformed version of a very well known viable inflationary model, namely the hilltop 
inflationary potential [93-100]. Thus, in a next subsection we will be able to provide a 
consistent framework consisting of two scalar fields for the aforementioned inflationary 
cosmologies. 


2.1 Singular evolution with hilltop-like potentials: Description with a single 
scalar and the singularity appearance 


We start our analysis with the hilltop potentials, and we demonstrate how a Type IV 
singularity can be consistently incorporated in the theoretical framework of these models 
by using a very general reconstruction scheme for scalar-tensor theories. Before getting 
into the main analysis, it is worthy to give a brief description of the scalar reconstruction 
method we shall use. For a detail presentation on this issue the reader is referred to [1, 102], 
We consider the following non-canonical scalar-tensor action, which describes a single 
scalar field: 

S = J d?xy/=gl^R- ^oj^d^d^c/)- V(<t>) + T matter | , (2.3) 


with the function <u(</>) being the kinetic function and V((j>) the scalar potential. The 
kinetic-term function appearing in (2.3) is irrelevant, since it can be consistently 

absorbed by appropriately redefining the scalar field <f> as 


V = 



(2.4) 


where it is assumed that w(</>) > 0. Hence, the kinetic term of the scalar field appearing in 
action (2.3) can be written as 


- oj^d^d^cj) = -d^ipd^tp . (2.5) 

On the other hand, if cj(</>) < 0 then the scalar field becomes ghost, and it can describe the 
phantom dark energy. In this case, instead of (2.4) one makes the redefinition 

d(f>yj-bj{4 >), ( 2 . 6 ) 

and thus, instead of (2.5), one acquires the expression 


- = d^ipd^ip . (2.7) 

When the scalar field is a ghost field the energy density corresponding to the classical 
theory becomes unbounded from below. At the quantum level the energy could be rendered 
bounded from below, however at the significant cost of the appearance of a negative norms 
[108, 109], 

In both canonical and phantom cases, in the case of FRW geometry the energy density 
and pressure of the scalar field write as 

P = + V (<1 >), P = ~V{4>). (2.8) 
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Having these at hand, we can use the usual Friedmann equations in the absence of the 
matter sector, in order to promptly express the scalar potential V (</>) and the kinetic term 
uj((j)) in terms of the Hubble rate and it’s first derivative as 

= -4 H , V(<f>) = 4 (3 H 2 + h) . (2.9) 

The scalar-reconstruction method is based on the assumption that the kinetic term u((f>) 
and the scalar potential V(0), can be written in terms of a single function f{(f>) in the 
following way: 

= -4m > v(<t>) = 4 [ 3 /(^) 2 + f'wi • ( 2 - iq ) 

hi hi 

Therefore, the equations (2.9) become 

4 = t, H = f(t). (2.11) 

Finally, by varying the action in terms of the scalar field we obtain its evolution equation 
as 

0 = ui(<j>)(f> + -u/+ 3+ V'(4>) , (2.12) 

which is also satished by the solution (2.11). 

In order to proceed, we need to make a choice for the scalar potential. In the case 
of a canonical scalar field, the hilltop inflation potential [93-100] has been proved to be 
consistent with the Planck data [21], and it is approximately given by 


V{ip) + , (2.13) 

with the ellipsis denoting higher-order terms which are negligible during inflation. When 
the parameters are chosen appropriately, and amongst others p = 2, the implications of 
this inflationary mode are in 95% agreement with the Planck data [21]. Hence, we prefer to 
consider a case close to this phenomenologically well-behaved case, and we assume p ~ 2—e, 
with e an infinitesimally small number |e| <C 1. 

Let us now start the main part of the analysis of this subsection, which is to examine 
whether there exist finite-time singularities in the hilltop potential scenario (2.13). As we 
shall demonstrate, it is possible to connect a Type IV and a Type I singularity with poten¬ 
tials having a functional resemblance to the hilltop potential. However, it is impossible to 
fully satisfy the constraints posed by Planck in these hilltop-like potentials, and therefore 
we should stress that the following considerations are excluded by observations. Neverthe¬ 
less, for the moment we are interested in demonstrating that instabilities will appear, and 
we do not consider consistent confrontation with observations. One could construct more 
complicated single-field scenarios, in agreement with observations, and with the appear¬ 
ance of singularities, however this would lie beyond our scope, which is to show in a simple 
single-field model that singularities appear. 

In order to quantify the appearance of a singularity we parameterize the Hubble rate 
as 

H(t) = co + b 0 (t s - t) a , (2.14) 
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where Co, &o are arbitrary positive constants, and later on we shall determine which values 
of a are allowed. Clearly, according to the value of a , one could have the appearance of 
the following singularities: 

• a < — 1 corresponds to Type I singularity. 

• — 1 < a < 0 corresponds to Type III singularity. 

• 0 < a < 1 corresponds to Type II singularity. 


a > 1 corresponds to Type IV singularity. 


Inserting the above Hubble ansatz into (2.10), we may directly find the scalar potential 
and the kinetic term as 


u{4>) = 


— 1+0! 


V{4>) =-$- - 


2 b 0 a(t s - (j)) 


K 

3eg b 0 a(t s - (j)) 


— 1 -\-a 


+ 


6b 0 c 0 (t s - 4>) a 3 bl{t s - 4>) 


2a 


+ 


K“ K* K 

Additionally, imposing the transformation (2.4) we obtain 


ip = 


2 ' /S “ (;t 


k (a + 1) 

and therefore the scalar potential (2.16) becomes 


V(<p) ~ - 


3cg b^aAi 


a— 1 


2 (- 1 + a ) 660C0A? 36 nA? Q 

ip o + l --- L(pa+1 -|- U 1 + l 


AC 


with 


Ai — 


2y/2oibo 


2 

0 + 1 


(2.15) 

(2.16) 

(2.17) 

(2.18) 

(2.19) 


k (a + 1) 

Since we are interested in small values of ip, which according to (2.11) corresponds to 
the case where t s — 4> —> 0, we may neglect the higher-order terms in the potential (2.18), 
and approximate it as 




2(°-l) 

o+l 


* 

K AC 

In the case p = 2 — e this potential coincides with (2.13) under the identification 


a = 


2 + e 


( 2 . 20 ) 


( 2 . 21 ) 


which is a very large positive number as e —>• 0, if e > 0. Hence, in this case the system 
develops a Type IV singularity. On the other hand, when p > 2 the system develops a 
Type I singularity, since in this case comparing (2.13) with (2.20) we obtain 


a = 


p + 2 

2 — p 1 

which is clearly negative for p > 2 (i.e. e < 0). 


( 2 . 22 ) 
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Before proceeding, let us make a comment on the above reconstruction procedure. 
Actually, this procedure is based on the Hamilton-Jacobi-like approach independently in¬ 
troduced in [110, 111], however now the difference is that instead of solving the resulting 
first-order nonlinear differential equation for H(<f>) (or equivalently H{t )) we impose it 
at will, and we calculate the potential V(4>) required for the given dynamical evolution. 
Strictly speaking, one should examine whether there are other solutions corresponding to 
the same potential, but for the purposes of this work it is adequate to extract one of them. 
We now calculate the effective equation of state 


P 1 

p 3 H 2 ' 


Using the Hubble parameter (2.14) we obtain 


(2.23) 


w e & = -1 + 


2 b 0 (-t + t s ) 1+a a 
3 [cq + bo(—t + t s ) a ] 


(2.24) 


In one of the interesting cases at hand, namely when q > 1, which corresponds to the 
Type IV singularity, the effective equation-of-state parameter becomes exactly equal to 
iCeff = — 1 at the time t = t s . This case is problematic for the following reason: As one 
can see, the main source of the problem originates from the fact that according to (2.14) 
we obtain H = 0 at t = t s . This can potentially introduce a very large instability when 
crossing the phantom divide w e g = — 1. In order to show this explicitly we introduce the 
variables and Y , defined as 


X,.p 



Y 

H 


(2.25) 


In practice, the variable Y quantifies the deviation from the solution of the reconstructed 
scalar-tensor theory, given in (2.11). Using these new variables, the Friedmann equations 
(2.9) and the field equation (2.12) can be rewritten as 


dX^ _ f'W(x%-l) dV _ f (</>)( 1-AVU)X 0 

d N 2 1 * ’ ’ dN H 2 


(2.26) 


where N is the e-folding number. Hence, the scalar-tensor reconstruction solution of (2.11) 
corresponds to the following values for the new variables: 


X 0 = 1, Y = l, 


(2.27) 


which correspond exactly to the basic critical point of system (2.26). As usual, in order 
to examine when this solution is stable we linearly perturb it around this critical point as 
[112-114] 

X^ = l + 5X (j) , Y=1 + 5Y, (2.28) 

and thus the dynamical system (2.26) can be re-written as 

J-&-* 

\ 




d 

d N 


(2.29) 








Hence, the stability of the dynamical system (2.29) is ensured if the eigenvalues of the 
involved 2x2 matrix are negative. These eigenvalues in general read as 


M + 


1 

2 


( H"{t) H'(t) \ 


+ 


f H"(t ) H'(t) \ 2 4 H"(t) 12 H'(t) 

\H'(t)H(t) + H (t) 2 + 3 J H(t) 3 H(t)' 2 




1 

2 


/ ff"(f) if'(t) 


f H"(t ) if'(t) \ 2 4H"(f) 

(fr(^) + W + 3 J 


12H'(f) 

~hW 


and thus for the Hubble rate (2.14), at the time t ~ t s , they become 


M+ = 0, 




2(—1 + a ) 
(~t + t s )c 0 


(2.30) 


(2.31) 


Since M_ > 0 when t < t s and M_ < 0 when t > t s , the dynamical system is unstable 
at the transition point t = t s , which corresponds to the phantom-divide transition point. 
This result is valid for any value of the parameter a. In addition, at the transition point 
the eigenvalue is infinite. This instability was first observed in [18, 115]. In order to 
evade this kind of problems in scalar-tensor cosmologies, the reconstruction method must 
be enriched with the presence of two scalars. As we show in the next subsection, the 
instability does not occur in such a case. 


2.2 Singular evolution with hilltop-like potentials: Description with two scalar 
fields 

As we demonstrated in the previous subsection, when we consider the transition from the 
non-phantom phase to the phantom one, in the context of scalar-tensor gravity with one 
scalar field, there appears an infinite instability at the transition point [1, 18, 102, 115]. 
However, this instability can easily be removed by considering scalar-tensor theories with 
two scalar fields. 

Consider the following two-field scalar-tensor gravity, with action: 

S = I d 4 x^g | ^R - - ^r/(x)<9 M x^ M X - V(<j>, x) j ■ (2-32) 

As in the single scalar case, the function uj(4>) is the kinetic function for the scalar field 
(f> , and accordingly rj(x) is the kinetic function for the scalar field x■ When the functions 
or rj(x) are negative the corresponding scalar field becomes a ghost field. However, 
since such a field would lead to inconsistencies at the quantum level [108], one expects it 
to arise through an effective description of a non-phantom fundamental theory [116]. 

Assuming that the two scalar fields 4> and x depend only on the cosmic-time coordinate 
t, and also that the spacetime is described by a flat FRW metric of the form (2.1), we easily 
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write the Friedmann equations in the form 


(2.33) 

(2.34) 


*m 2 +v(x)x 2 = , 

v^x) = ^(zH 2 + h) . 

Following the procedure of the previous subsection, and generalizing it in the case of two 
fields, we impose the parametrization 

w(f) + rj(t) = -^2 fit) , V(t) = \ [3/(f) 2 + /'(*)] , (2.35) 

K Kj 

which is consistent with an explicit solution for (2.33) and (2.34) of the from 

<t> = X = t, H = (2.36) 


For a detailed account on this reconstruction method see [1]. 

The kinetic functions ui((f>) and r/(x) can be chosen in such a way that is always 
positive and rj(x) is always negative, hence one field is canonical and the other one is ghost. 
This is exactly the realization of the quintom scenario [4]. In particular, a convenient choice 
for the kinetic functions is 

v(<f>) = - ^2 {/'(0) “ V a i(4>) 2 + /'WO 2 } > 0, 

^(x)=-4v / «i(x) 2 + /'(x) 2 <0, (2.37) 


where a\{x) is an arbitrary function of its argument. We now define a new function /(<(>, x) 
through 


/(</, x) 



J dct>uj((f>) + J 


dxv(x) , 


(2.38) 


which has the characteristic property 


f(t) =/(*)» 


(2.39) 


and therefore the constants of integration are fixed in (2.38). Assuming that V((f>,x ) is 
given in terms of the function /( cj>, x) as 


mx) 


i 

^2 


3/(</>, X) 2 + 


#/(<£. x) 

d<f> 


x) 

dx 


(2.40) 


we find that, in addition to the Friedmann equations given in (2.33) and (2.34), the following 
two-scalar held equations are also satisfied: 

0 = u(4>)<j) + + 3 Hu((j>)(t> + ^ ; (2-41) 

0 = n(x)x + \v'{x)x 2 + 3 Hr](x)x + ^ • (2-42) 

2 dx 
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In summary, upon using the above kinetic functions uj((f>), r/(x) and the scalar potential 
V { we have a two-scalar field scalar-tensor model with cosmological evolution of the 
form given in (2.36). 

Let us now extend the analysis of the previous subsection, and proceed to the recon¬ 
struction of a two-field scalar-tensor gravity generating the cosmology described by the 
scale factor (2.14). As we discussed, this class of cosmological evolution is related to a 
single-field scalar-tensor theory, which generates a hilltop-like scalar potential of the form 
(2.13) for the corresponding canonical scalar held. Moreover, it exhibits specific finite-time 
singularities, according to the values of a given in the list below (2.14). 

In order to proceed, we need to impose an ansatz for the arbitrary function oi\{x) 
appeared in (2.37). We choose it as 

ol\(x) = \/b 0 a(-x + tg) 01 - 1 - H'{x) 2 , (2.43) 


and therefore (2.37) give 


u{4>) = - 

v(x) = ~ 


b 0 a(t s - (f>) 1+01 + v 7 boa(t s - (f>) 1+c 


2^b 0 a{t s - x)~ 1+a 


At 


In order to avoid inconsistencies we assume that a has the form 

n 

a = 


2m + 1 


(2.44) 


(2.45) 


and therefore when n is an odd integer, the kinetic function uj(cj)) is positive, while ij{x) is 
negative. The corresponding function /(</>, x) from (2.38) writes as 


f{<t>,x) = ~bo(ts ~ <t>Y - 


2 x%a(t s - 4 >) 1 2 \fboa{t s - x) 1 


{ex + 1) 

and thus the scalar potential V{(f>,x) hi (2.40) becomes 


(« + 1 ) 


(2.46) 


T/Y , X b 0 a{t s -(j)) a 1 2^%a{t s - (j))° , 2 + x%a{a - 1) {t s - <j>) “2 

V{(p,x) = - 5 - 1 - 


(1 + a) At 2 


+ 


^_ Q, _ ^ ^_ 

2 Vboa(t s -x)— >/boa{a - l)(i s - x) 


Oi — 1 


(1 + a)t c 2 


(1 + o)k 2 


+ 


b 0 (t s - <A)“ + + 2V ^ {t - - X) ‘ 1 


(a + 1) 


(a + 1) 


(2.47) 


For the above two-held scalar-tensor theory, the solution (2.36) is stable as we will 
demonstrate in detail. Indeed, we introduce the quantities X^, X x and Y through 




*x = X, Y = 


!{</>, x) 

H 


(2.48) 
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and thus the Friedmann equations (2.33), combined with (2.41) and (2.42), can be expressed 
as 


dx^ _ (x* - l) 

dN 2u((j))H 

dX x _ rf{ X ) {X 2 x - 1) 
dlV 2 r,( X )H 

dV _ 3X^X x (1 - Y 2 ) 

div x (P + x x 


-3 {Xf-Y), 

-3 (X x -Y) , 

H X^X x + 1 - Y(X ep + X x ) 
H 2 V 0 + 


(2.49) 


The solution appearing in (2.36) corresponds to the following values of the variables X^, 
X x and Y: 

X 4> = 1, X x = l, Y = 1, (2.50) 

which correspond exactly to the basic critical point of system (2.49). As we did in the 
previous subsection, in order to examine when this solution is stable we linearly perturb it 
around this critical point as 


= 1 + 6 X 4 , X x = l + SX x , Y = l + SY, 


(2.51) 


obtaining the perturbative dynamical system 


d 

dN 


( SX A 

5X x 

V ft ) 


( 

V 


_ o 

Hw{<!>) 0 

0 

0 


0 

v'ix) 

Hv(x) 

0 


3 ^ 

3 

o H 

J 



(2.52) 


The eigenvalues of the matrix appearing in the dynamical system (2.52) are found to be 


w '(0) o v'(x) o A/r o H co , 

Hu :(</») 3 ’ Mx ~ Hn(x) ’ Y ~ # 2 ' (2 ' 53) 

Hence, in the case where the Hubble rate is given by (2.2), these eigenvalues become 

(-1 + a) (l + 2 y/b 0 a(t s - <^)~ 1+Q! ) ^b Q a(t s - (j))- 1+a 


Mih = —3 — 


2co 

— 1 Qt 


b 0 a(t s - (j)) a + (t s - 4 >)yjb 0 a{t s - 4 >) 1+01 


M x = - 3 - 

2 c 0 {t s - x) 

My = -3 , 


(2.54) 


which are clearly negative for all values of the cosmic time t. Therefore, we deduce that the 
solution (2.36), with the Hubble rate given in (2.2), is stable. Hence, as we have already 
mentioned, the presence of two scalar fields is adequate to solve the instability problem 
of the one-field theory. Having the stable solution (2.36) at hand, we establish the fact 
that a Type IV singularity can underlie the cosmology generated by the Hubble rate (2.2). 
In particular, one possible description can be realized by a two-field generic scalar-tensor 


- 12 - 






















theory with kinetic functions given in (2.44) and with the scalar potential appearing in 
(2.47). 

Finally, it is worthy to mention that the solution (2.36) is only one particular solution 
of the dynamical system (2.33) with (2.41) and (2.42). In principle, there might exist other 
solutions, however the stability of (2.36) which renders it an attractor, ensures that this 
solution will be the asymptotic limit of a large class of solutions. 

2.3 Singular versus non-singular evolution 

Let us close this section by making some comments on the possible removal of the singu¬ 
larities discussed above. In the literature there are mechanisms that could remove some 
singularities of a given cosmological scenario. Two of the known such mechanisms are the 
inclusion of quantum effects, or the suitable slight modification of the action (for instance 
changing the F(R ) or the potential form [117]). 

Indeed, singularities can become milder through the incorporation of quantum effects 
in a perturbative way, without changing the action. However this is not the case for the 
singularities analyzed in the previous subsections (in general even the quantum description 
of the phantom regime is ambiguous [108, 109]). 

Nevertheless, one could regularize and remove some of these singularities by modifying 
the action, i.e. by slightly changing the potential [117]. However, following this approach 
is not the scope of our analysis. In particular, we would like to mention that in principle 
the potential of a model is expected in general to be determined by the fundamental theory 
itself, thus in general one cannot modify it at will in a suitable way that will make the 
singularities to be absent. Hence, since any potential could arise from a fundamental theory, 
even if it is of special form, including the ones considered in this work, it is interesting and 
necessary to investigate their cosmological implications and in particular the realization of 
singularities, without entering the discussion to modify it, i.e modifying the theory. Indeed, 
a theory admitting weak singularities has to be special in some respect (for instance note 
that there are not such singularities in the vacuum general relativity case), however this 
could still be the case in Nature, deserving further investigation. 

Additionally, note that there could still be cases where weak singularities cannot be 
removed even after a slight modification of the action. For instance, as it was shown in 
[118], the singularities arising in F(R) gravity in the case where F"(Rq) = 0 for some 
R = Rq, may not be removed through the procedure of [117] if F"'(Rq) / 0. 

In summary, although removing the singularities is definitely an option, and a non¬ 
singular evolution can be realized, since singular evolution is a possibility too in this work 
we are interested to investigate the resulting cosmological behavior. Moreover, we desire 
to show that in some cases a singular evolution can be consistent and in agreement with 
the standard observed Universe history. 

3 Two-field scalar-tensor gravity: specific examples 

The use of two scalar fields offers a great number of possibilities for cosmological evolution, 
with various diverse features that were absent in the case of a single scalar field. In this 
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section we are interested in reconstructing specific two-field models that may still exhibit 
some singularity types, in particular Type II and Type IV. 


3.1 Generic examples of two-field scalar-tensor theories exhibiting Type II or 
Type IV singularities 

Let us start by considering the Hubble rate to have the very general form 

H(t) = h(t) + f 2 (t)(t s -t) a . (3.1) 


We can choose the arbitrary functions fi(t) and / 2 (f) to be smooth and differentiable 
functions of t. In the following, and without loss of generality, we shall consider the case 
where a is given by two integers n and m, namely 


n 

2m + 1 


(3.2) 


and specify furthermore the value of the integer n when this is needed. For more general 
values of the parameter a we may extend (3.1) to 


H(t) = fi(t) + f 2 (t) |t s - t\ a . 


(3.3) 


Concerning the kinetic functions to(4>) and r)(4>) of the two fields that appear in action 
(2.3), we use 


w(<£) = » v(x) = |/2(x) (t B ~ x)° + af 2 (x) (t s ~ XT 1 

Kj t\j 

which implies that the function a\{x) appearing in (2.37) is of the form 


ai (x) = \l f 2 (x) (t s - x) a + af 2 (x) (t s - x) a 1 ~[H'(x)} 


i 2 


(3.4) 


(3.5) 


Here, the variable x can be either the scalar field cj) or X■ Having chosen ot\{x) as in 
expression (3.5) allows to specify the final form of the function f(<j>,x), defined in (2.38), 
and indeed in the case at hand becomes 


f(<t>,x) = /i(<£) + f“i (x) (t s - x) a ■ 


(3.6) 


The corresponding scalar potential of (2.40) reads 


mx) = ^2 { 3 + h(x) (ts - xY } 2 + f[( 4 >) + Y(x) (t s - xY + a fz(x) (t s 

K\ v 

Additionally, the kinetic functions from (3.4) now become 


x)”- 1 } . 

(3.7) 


u(4>) = 

v(x) = 


2 m 


2 [a(t s - x) 1 +a f2(x) + (ts ~ xYY(t )] 


(3.8) 
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Finally, from the form of the Hubble rate (3.3), we can straightforwardly extract the total 
equation-of-state parameter of the universe (2.23) as 


w e s = -1 - 


2 + (~t + t s ) a f 2 {t) ~ (~t + t s ) 1+a af 2 (t)] 

3[/i(i) + (~t + t s ) a f 2 (t)\ 2 


(3.9) 


As we already mentioned, the functions fi(t) and f 2 {t) are arbitrary functions, so a 
convenient choice might lead to interesting results, regarding the cosmological evolution. 
Assuming that the Hubble rate is equal to H(t ) = fi(t) + f 2 (t), one convenient choice for 
/i(t) and f 2 (t) could be 

hit) = f^), f 2 (t) = f x (x), (3.10) 

where 


fHd>) = f}(<t>) + fttt) (h ~ </>)“ , f x ix) = fhx) + fHx) (t 2 - xf , (3.11) 

and where /(f (</>), f\(x)> an d f 2 {x) are smooth functions of the scalar fields. Then, 

in terms of the functions appearing in (3.11), the Hubble rate is written as 

H{t) = ft’it) + f x (t) = ff(t) + ff(t) + f$(t) (ti - (/>)" + / 2 x (t) (t 2 - . (3.12) 


An interesting cosmological scenario can be realized if we choose t\ <C t 2 , a > 1 and 
(3 < 0. In this case, t\ could correspond to a time around the inflation era, while t 2 to the 
dark energy epoch, present or future. In addition, when a > 1 and /? < 0 it easily follows 
from the Hubble rate (3.12) that at t\ the physical system of the two scalars experiences a 
Type IV singularity, while at t 2 , it experiences a Type I (Big Rip) singularity, if j3 < — 1, or 
a Type III, if —1 < (3 < 0. It is worthy to further specify the functions appearing in (3.11) 
and quantitatively examine the cosmological evolution and the resulting physical picture. 
As an example we choose 

f*W = T , + Tpfn (*i - *)° , f x (x) = h(t 2 -xf , (3.13) 

V t 2 0 + </> 2 tg + 

where to, f\, f 2 , and fs are considered to be positive constants, so that H(t) > 0. For the 
choice (3.13), the Hubble rate becomes 

and the corresponding kinetic functions of the two scalars read 


ui4>) 

vix) 


K 


fit 


(ti + 4f 2 
2 hi-t + t 2 )~ l+ ^ 


4f 2 t 5 (-t + h) a 2f 2 t(-t + h) a 

{t 4 + 4) 2 t 4 + t* 


f 2 t 2 (-t + ti) 1+ol a 

t 4 + 4 


(3.15) 
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Since the parameter f3 is assumed to be /3 < 0, the kinetic function rj(x) is negative, and 
therefore x is a ghost field. On the contrary, u(cf>) is positive and therefore <f> is a canonical 
scalar field, as long as / 2 -C f\. The corresponding f(cj>,x) function from (3.6) writes as 


fvt>, x) = ,... + — tt—ta — + h\P 2 - xr 


t 0 + ■ 


\Ao + <\> 2 

and therefore the scalar potential (3.7) is given by 


V{(j),x) =\ 


fi<t> 


(^ + 0 2 ) 3/2 (4 + ^) 


4/ 2 (ti - 0)“^ 5 2 Mh - </>)*</> / 2 a(ti - </)- 1+Q </ 2 


tn + <^ 4 


tn + 


+ 


/1 , -/>)“</ f ^ 
7fr? + +/3<2 ' x) 


- - x) 1+/? 


(3.17) 


Finally, the total equation-of-state parameter corresponding to the Hubble rate (3.14) reads 


2 

w eS = -1- 


_ fit 

{t 2 +tl) 


Tj2 


4/2i 5 (-Wi) a , 2fjt{—t+t\) a _ fit 2 (—t+ti) 1+a g 

(t 4 +t^Y H+t 0 t +t 0 


3 


A , f2t 2 (-t+t 1 y 


2 


+ /s(-^ + ^) /3 


f3(~t + t 2 ) 1+/3 /5 


(3.18) 


In the early universe /^(t) dominates in H(t ) in (3.12), and therefore the inflationary 
era occurs, with a Type IV singularity at t = t\. In the late-time universe f x {t ) dominates 
in H(t) in (3.12), and the universe will end with a Type I (Big Rip) singularity (j3 < — 1) 
or Type III singularity (—1 < f3 < 0) at t = t 2 . Having made these assumptions, it is 
worthy to study the behavior of the EoS (3.18) as a function of the cosmic time. Let us 
first consider time scales near the Type IV singularity, thus t ~ t\. In such a case, by 
disregarding subdominant terms, the EoS takes the form 


rc eff ~ -1 - 


- 7 ^ 572 -/ 3 ^ 2 ) 1+/3 /3 

(*o) 


% + h<ti) p 


(3.19) 


and since t 2 1 and /3 < 0 it finally becomes 


ICeff ~ -1 + 


2 1 

3/Ho ’ 


(3.20) 


which describes quintessential inflation. In the case where to 1, the EoS (3.20) becomes 
approximately w e s — —1, i.e. it is effectively a cosmological constant. Furthermore, with 
regards to the late-time behavior, which corresponds to cosmic times near t 2 , we shall 
assume that t 2 > to and also that a is of the form given in (3.2), namely a = n/{2m + 1). 
For n even, and 2 < a < 3, the effective equation of state reads 


t , 2/ 3 (-t + t 2 )- 1+ ^ 

W e g ~ — 1 -i-—s , 

3lf2t~ 2+a + f3(-t + t 2 )P} 2 


(3.21) 
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which is Wes < — 1 when t < t 2 . Therefore, eventually this case describes phantom dark 
energy at a time before the singularity. Additionally, the case a > 3 and n even is slightly 
more complicated. In this case we set x = — t + t 2 , and the EoS becomes 


ICeff — -1 - 


2(t 2 — x) 1 ~ a {—2 + a) 
3/2 


(3.22) 


which for t < t 2 (or x > 0) is w e s < — 1, which again describes phantom evolution. In the 
case where n is an even integer and a < 3, the EoS reads 


, , 2/ 3 (—f + t 2 )- 1+ ^ 

W e tt — —1 H-n , 

3[f2t- 2+a -M-t + t 2 )P} 2 


(3.23) 


which describes phantom acceleration (recall that f3 < 0), while for a > 3 we have 


w eS ~ -1 - 


2(* 2 -x) 1 ~ a (2 + a) 

3/ 2 


(3.24) 


which describes phantom acceleration too. In summary, with two scalar fields we may 
realize a plethora of cosmological expansions, with finite-time singularities present in the 
cosmological evolution. 

We close this subsection with another model with an interesting cosmological evolution. 
In particular, we consider the Hubble rate to be 


H(t) = ff(t) + /*(£) + f 2 (t) (t s - t) a . (3.25) 

The function ff(t) can be chosen in such a way that it describes the (smooth) inflation era, 
and in addition the function /*(t) can be chosen to describe the present time accelerating 
expansion of the Universe. We can consider the following two classes of models, in which 
the functions are: 


(A) : /^) = fttf) + M</>) {t 8 - 4>) a , /*(*) = f?(x ), (3-26) 

(B) : /*{</>) = , /*(x) = f?(x) + MX) {ts ~ XT • (3-27) 

Obviously, both models (A) and (B) generate identical evolution of the Universe’s expan¬ 
sion, since the Hubble rate is the same. In model (A), the scalar field 4> can be selected in 
order to generate both the inflationary era and the finite time singularity. In addition, in 
model (B) x can generate both the present accelerating expansion, as well as the finite time 
singularity. We may choose t s to be large enough, in order to correspond to the far future, 
and moreover the function f 2 (t) can be chosen in order to give a small contribution to the 
present cosmological evolution of the Universe, and hence it is difficult for it’s contribution 
to be observed. In the following, we shall present a model for which this scenario can be 
realized. 

For brevity, let us consider model (3.26). The Hubble rate is simplified as 

H{t) = f?(t) + {-t + h) a f$(t) + fi(t ), (3.28) 
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and the kinetic functions of the scalar-tensor theory read 


u>0) = -- 


-i-t++rt(t) + (~t+h 


while 


v(x) = ~ 


K 2 

-± , (3.29) 

2 m) 

k-2 ' 

(3.30) 


The corresponding scalar potential (3.7) becomes 

=- I (*i ~ + 3 \ft{4>) + (h ~ + fiix) 


while the total EoS writes as 


u>efi = -1 - 


-(-t +ti)- i+ “«/ 2 ^(t)+ ft® + (~t + hr?*® + mt) 


3 [flit) + (-t + h) a f$(t) + /f(t) 

A convenient choice for the functions appearing in the Hubble rate (3.28) is 

ft(t) = —JL=, f+(t)=e-v, m- ht2 


(3.31) 


(3.32) 


V t2 + t o 

and in this case we finally acquire 
.fit _ 4 h£ 


t 4 + td’ 


(3.33) 


u>es = -i - 


(t2 +i gp (t*+tg) 


^ + t/3 (-£ + U) 1+Q a-e t/3 (-t + ti) Q /3 




/2* 2 


(3.34) 

Thus, we can now study it’s behavior, focusing mainly at times corresponding to the 
inflationary era, as well as at late times. 

We first consider that t\ corresponds to the inflationary regime, and that at this point 
a Type IV singularity occurs. This implies that a > 1 and moreover we assume that (3 > 0. 
Furthermore, if to > U then the EoS at early times reads 


w eS = -1 + 


2f(—2/ 2 + /it 0 ) 


2+2 


3,/ft, 


(3.35) 


In the case where /Ho > 2/ 2 , the Universe accelerates in a quintessential way near the Type 
IV singularity at t\. In addition, if to 2> 1 then the second term in (3.35) is negligible, and 
therefore w e g ~ — 1, which corresponds to a de Sitter expansion. Definitely, if /Ho < 2/ 2 
then the EoS (3.35) describes phantom inflation. 

Let us now focus at late times, where the effective EoS behaves as 


rc eff ~ -1 + 


21 


3/i\A 2 + t o ’ 


(3.36) 
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which describes quintessential acceleration. In the case where t ~ to, and to 1, relation 
(3.36) becomes a nearly de Sitter late-time acceleration, with w e g — —1. It is worthy to 
say that the scalar field cj) is the one with the positive kinetic term in all cases, with 
being equal to 


u{4>) = 


2 


( 2 f %/2 ~ e + ti) 1+a a — e ^(-t + ti) a P 

(t 

k 2 


(3.37) 


while r](x) reads 


4/2 ( t 5 - ttp 

{t 4 + 4) 2 K 2 ’ 


(3.38) 


and since to 1 it can be negative for a wide range of t values (we assume that f<i > 
0). Therefore, the use of multiple scalar fields can significantly increase the cosmological 
scenarios for the physical system described by a specific Hubble rate. 


3.2 Slow-roll with two-field scalar-tensor theories 

In this subsection we qualitatively describe how the slow-roll approximation is affected by 
the use of two scalar fields describing the cosmological evolution. In addition, we discuss 
how the singularities affect the general slow-roll evolution and also we investigate when 
these can lead to inconsistencies. In principle, the results can be model-dependent, but 
using some illustrative examples we shall see how a singularity can affect the inflation 
parameters. Our study will remain at a qualitative level, since a detailed investigation 
would require the analysis of two- and three- point statistics of field perturbations, which 
can be related to the two- and three- point statistics of the curvature perturbations on 
uniform density hypersurfaces. This could be achieved by using the <5N-formalism [119— 
127], however this detailed study lies beyond the illustrative scopes of this article. For 
details the reader is referred to [119-127]. 

The most important inconsistency that the singularities might cause, is traced in the 
definition of the e-folding number N. The inconsistency occurs when the singularity ap¬ 
pears during the inflationary era, even in the case the singularity is of Type IV. Therefore, 
we shall assume that the singularity occurs at exactly the time when inflation ends. In the 
end of this subsection we shall discuss in detail this problematic issue. 

In general, single-field and multi-field inflation are very different, when the correspond¬ 
ing scenarios are viewed as dynamical systems. In the single-field case, the phase space 
formed by the slow-roll solution is one-dimensional, implying that the slow-roll solution is 
the unique attractor of the dynamical system. On the other hand, in the case where two- 
field inflation is studied, the phase space is two dimensional and the stable trajectories form 
an infinite countable space [119-127]. Hence, the field values at the end of inflation will 
generally depend on the choice of the classical Held trajectory. The corresponding e-folding 
number will determine the evolution along one of these classical trajectories [119-127]. In 
order to see this in a qualitative way, in this subsection we shall investigate when the slow- 
roll condition is violated in the cases of single-field and two-field inflation. The violation of 
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the slow-roll condition will actually determine when inflation ends. We quantify our results 
by studying the behavior of the slow-roll inflation parameter e. 

As a toy single-field example we shall consider a viable inflationary model, which be¬ 
longs to the class of large-field inflationary realizations. In a previous article we performed 
a thorough analysis on this type of scenarios [20] , and thus we now focus on the qualitative 
analysis of the slow-roll parameter e. The model consists of a canonical scalar field p, and 
the potential for large held values is given by 

V(cp) ~ d(p - ti) n , (3.39) 

with t\ the time at which the Type IV singularity occurs, which for consistency is chosen 
to coincide with the time that inflation ends. As was demonstrated in [20] the Type IV 
singularity can be consistently accommodated in the cosmological evolution of the model 
under study. The slow-roll parameter e in the case of a single canonical scalar held reads 
as 

1 

e ~2^ 

and thus for the potential (3.39) it becomes 

Ti^ 

^ = 2 K?(t s -pY 

It is easy to analytically calculate when inflation ends, which happens at a held value p e for 
which e(p e ) ~ 1. Indeed, the slow-roll condition is violated when the scalar held exceeds 
the critical value <p e , namely 

Ti 

ip e — 1\ = • (3.42) 

V2k 

In order to acquire an overall picture of the phenomenon, in Fig. 1 we depict the behavior 
of e(tp) as a function of p. Notice that the large held values correspond to early times, and 
as the held value decreases, the cosmological time increases. Hence, the x-axis should be 
viewed backwards, as time increases. As we can see, the slow-roll condition is violated at 
the held value p e = t\ + which is indicated by a dot on the graph. Thus, in the single¬ 
scalar case the ending of slow roll is characterized by a single point in the corresponding 
phase space. This fact has a direct impact on the slow-roll parameter e(p), as can be seen 
in Fig. 1. 

Let us now come to the two-held case, which as we show is more involved. In order to 
simplify things we choose a simple example, and we consider the Hubble rate to be 

H(t) = fi(t) + f 2 (t), 

with fi(t) and / 2 (f) being equal to 

fi(t) = vi (~t + h) a , f2(t) = V2(-t + t 1 ) a , (3.44) 

where u\ and ^2 are parameters. In this case the Type IV singularity occurs if a > 1. 
Thus, we can accommodate the Type IV singularity in a generalized two-held scalar-tensor 


(3.43) 


(3.41) 


V(p) 


(3.40) 
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Figure 1. The slow-roll parameter e(tp) as a function of the scalar field ip, for the potential (3.39) 
for n ~ 1.95. As the scalar field decreases the cosmological time is assumed to increase. The dot 
corresponds to the critical field value <p e — t\ + , at which the slow-roll condition is violated. 


theory with action given in (2.32), by following the reconstruction method of the previous 
section, with a\{x) chosen as in (3.5). In particular, the kinetic functions oj(f>) and rj(x) 
are equal to 

w(<£) = —%/i( 0 ), v(x) = --^/2(x), (3-45) 

and by substituting the functions from (3.44) we acquire 


w(</>) = 


2 v\(—t + ti) 


a—l 


a 


v(x) = 


2v 2 {-t + ti) a l a 


(3.46) 


Therefore, the corresponding scalar potential V((f), x) writes as 

Tr/J . ^ viafa - <f>)~ 1+a 3 vf(ti-(j)) 2a v 2 a{t\ - x)~ 1+a 
V (0, X) =- 5 -1-o-o- 


K, 


K, 


GviMti - 0)“(*1 - x) Q 3 z/|(fi — x) 2a 

I Q I o 




(3.47) 


We assume that t\ is identical with the cosmological time when inflation ends, and thus it 
is a very small number. In addition, we assume that inflation occurs for large field values, 
and thus at the time when inflation occurs, t\ — cj) 1 and t\ — x 1- Since a > 1, the 
potential can be simplified to a great extent by keeping the dominant terms, becoming 


N = 3iq 2 (fi - f)) 2a + 6i/i^a(ti - (j)) a {ti - x) a + 3t/|(ti - x) 2a 

k-2 l-2 


(3.48) 


As a next step, we can transform the above two-held scalar-tensor theory to one that 
contains only canonical scalar fields. Indeed, by making the transformation 


t\ — a = 


rx _ 

/ \AKxRy ; 


h-<p= / y/u((t>)d(t), 


(3.49) 
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we can extract the identifications 


2 2 
- 4> + h = ai (ip — t-i) “+i , - X + h = «2 (cr — ii) “+ 1 


(3.50) 


where aq and «2 stand for 



2\J2av\ 

CX+1 


2\f2otV2 

aq = 

(1 + a)n 

5 

OL2 — 

(1 + a)n 


2 

ct+1 


(3.51) 


Hence, in terms of the canonical scalar fields <p and a, the simplified scalar potential (3.48) 
reads 


O O rv ZOL ZOL 

3vf(ai) 2a (ip - H)“ +1 , 61/1^2 (ctia 2 ) (<£ - H) a+1 (<? - H) “+ 1 
V (p, a) = -o-1- n - 


K 

3z/|(o;2) 2a (cT - ti)^+i 


+ 


(3.52) 


while the canonical two-field action is simply 


■ f 1 


1 


S = j <rxy/=g | - -d^ipd^ip - -d^ad^a - V(tp, a) 


(3.53) 


Having the scalar potential (3.52) at hand, we can qualitatively investigate when the 
slow-roll approximation is violated. Practically this occurs when the slow-roll parameter e 
becomes e > 1. In the case of a scalar-tensor theory with two canonical scalar fields, the 
slow-roll parameter e is equal to [119-127]: 


e(tp, a) = + e a , 


(3.54) 


with and e CT reading as 


1 

'd<pV(tp,ay 

2 1 

d„V (ip,cr)' 

^~2^ 

b 

_1 

’ e °~2^ 

V {ip, a) _ 


(3.55) 


Thus, by using the potential (3.52) we may directly evaluate the slow-roll parameter e 
and study it’s evolution as a function of the two scalar fields. Assuming again large field 
values for small cosmological times, which decrease as the cosmic time approaches the end 
of inflation, in Fig. 2 we present the behavior of the slow-roll parameter e as a function of 
the scalar fields. As we observe the slow-roll approximation is violated for a combination 
of values of the two fields, as expected. 

Before closing this subsection, we shall investigate the impact of the singularities on 
the cosmological parameters of the two-field model. In particular, we are interested in the 
cases where the singularities may lead to inconsistencies. As we mentioned earlier, the 
main inconsistency is traced in the definition of the e-fold number N, which is defined as 


N = — 



(3.56) 
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Figure 2. The slow-roll parameter e(<p,cr) as a function of the scalar fields (p and a, for the choice 
(3-44) f or = v 2 and a ~ 1.1. As the scalar field values decrease, the cosmological time is assumed 
to increase. 

where i* is the Hubble exit during inflation. The final time tf corresponds to a uniform- 
density hypersurface, at which inflation is considered to end. By taking into account (3.43) 
and (3.44), the e—folding number writes as 

N = (Vi + V 2 ) jfoi ~ t f ) a+1 - {-t* + H) a+1 ] 

ot -\~ 1 

Hence, it is obvious that when a < — 1 the above expression is always singular when tf = t \, 
due to the last term. However, if the singularity is of Type IV there is no inconsistency, 
which is the main point of this subsection. 

Finally, as we mentioned earlier too, we note that the complete analysis of the e-folding 
number for the two-field case can be quite complicated, and a detailed investigation would 
require the use of the <5N approximation. However, such a study stretches beyond the scope 
of the present work, and the reader is referred to [119-127] and references therein. 

3.3 Singular inflation with singular dark energy 

In this subsection we are interested in constructing a cosmological scenario where singular 
inflation gives rise to radiation and matter eras, and eventually to a dark energy epoch 
that results in a Big-Rip singularity. In order to succeed this, and following the previous 
subsections, we use two scalar fields, and thus this is another indication of the capabilities 
of such a model. 


(3.57) 
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Having in mind the discussion on the singularity types of the Hubble rate (2.14), we 
deduce that a Hubble rate that captures the above features, namely a singular inflation 
with a type IV singularity and a dark energy epoch resulting to a Big Rip, writes as 

H ( t ) = c 0 + b 0 ( t s - t) a + d 0 (t B - t) J , (3.58) 


where a > 1 and 7 < — 1, with co, bo, do constants, and where t s and t B are respectively 
the time of the inflationary type IV singularity and of the Big Rip. We need to note that a 
is assumed to take values of (3.2), therefore complex values for the Hubble rate are avoided 
for the physically interesting interval t s < t < t B . Practically, since the denominator of the 
fraction appearing in (3.2) is an odd integer always, then this procedure is always valid. 
Definitely, there are extra complex conjugate branches too, but we only focus on the real 
negative one. 

Assuming as usual that the two scalar fields (j) and x depend only on the cosmic-time 
coordinate t, and also that the spacetime is described by a flat FRW metric of the form 
(2.1), we extract the Friedmann equations in the form 

a#)</> 2 + v(x)X 2 = , (3.59) 

V(ct>,x) = ^(^H 2 + H) . (3.60) 

Following the procedure of the previous subsections we impose the parametrization 

w(i) +v(t) = > 

V(t) = ±[3f(t) 2 + f\t)] , (3.61) 

which is consistent with the explicit solution of (3.59) and (3.60): 

<l> = X = t, H = (3.62) 


A convenient choice for the kinetic functions would be 

w(<£) = --^2 {f'ifi) ~ V«iM r -KTWO 2 } > 0, 

V(x) = --^V«i(x) 2 + fix ) 2 < 0, (3.63) 


with <37 (x) an arbitrary function of its argument. It proves useful to define a new function 
f( 4 >, X) through 


x) 



J d(f>uj((f)) + J 


dxd(x) , 


(3.64) 


which has the characteristic property f(t) = f(t ) which fixes the constants of integration 
in (3.64). Assuming that V((j),x) is given in terms of the function f(cj),x ) as 




m , x f + ?MA + d J^A 


(3.65) 
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we find that in addition to the Friedmann equations given in (3.59) and (3.60) we have 


0 = u}((f>)4> + X ^ , (3.66) 

2 o<p 

o = v(x)x + \v'ix)x 2 + 3 Hr](x)x + ^ . (3.67) 

2 &x 

Hence, using the above kinetic functions co((f>), r](x) and the scalar potential V(cf>,x) we 
obtain a two-field scalar-tensor model with cosmological evolution of the form (3.58). 

As a next step we need to impose an ansatz for the function a\{x) appeared in (3.63). 
We choose it as 

«i(x) = \/b 0 a(—x + t s ) a_1 — H'(x) 2 , (3.68) 

and therefore (3.63) give 


u(4>) 

v(x) 


2 


b 0 a(t s 


(t>) a 1 + do'yits - cj)) 1 1 + y 7 bo<x{t s ~ <t>) a 1 

k 2 


2 ^/b 0 a(t s - x) 1+a 

K 2 


(3.69) 


As expected, one of the fields is canonical while the other one is ghost. The corresponding 
function f(<f>,x) from (3.64) becomes 


/(</>, x) 


-b 0 (t s - (j)) a - d 0 (t B - 0) 7 


2 y/b^a(t s - i))^ 1 

(« + 1 ) 


2 y/b^a(t s - x) a * 
(a + 1) 


(3.70) 


and thus the scalar potential V(4>,x) from (3.65) becomes 




b 0 a(t s ~4>) a 1 + d 0 x(t B - (t>V 1 

K, 2 

2 Vkja(t s - (j))^ + y/b^a{a - 1) (t s - 0)^ 
(1 + a)n 2 

2 Vb^a(t s - x)^ _ y/bpaja ~ 1 )(t s - xY -1 

(1 + a)n 2 (1 + a)n 2 

+ -4 { fe o(fs - (t>Y + do(t B - 

K 2 I (a + 1) 


. 2 Yb^(t s -x)^V 
(«+1) j ■ 

(3.71) 


In summary, the above two-field model can give rise to the Hubble rate (3.58), which 
corresponds to a universe with a singular (of type IV) inflation, resulting in a dark-energy 
phase that ends with a Big Rip. 


3.4 Barrow’s two-field model 

In this subsection, for completeness, we present the realization of the singular inflation 
in the Barrow’s two-field model [74], This model can exhibit the Type IV singularity 
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consistently incorporated in it’s scalar-tensor theoretical framework. In particular, we 
consider a Hubble rate of the form 


H{t) = b 0 it s - t) a . 


We choose the arbitrary function ai(x) appearing in (2.37) to be 

«i(x) = v 7 2otboi~x + t s ) (Q - 1)/2 , 


(3.72) 


(3.73) 


and substituting this into (2.37) we can easily obtain the corresponding kinetic functions 
uii<p) and rjix) as 


w(0) = 

v(x) = 


-b 0 ait s - 4>) 1+a - y/boa(ts 


hi 


2^b^ait s - <f>) * 2+ “ 



(3.74) 


In the case at hand, and assuming that a has the form (2.45), namely a = n/{2m + 1), 
with n constrained to be an even integer and a > 1, the kinetic function uj((p) is rendered 
positive, while the kinetic function rjix) becomes always negative. Thus, we can promptly 
find the function fi<p, x), which in this case is 




1+a 


-boil + a){t s - <P) a ~ 2y/%ait s - cp) 2 


(1 + a)n 2 


+ 


4y/%ait s - x) 1+2 ‘ 1 


(1 + a)K 2 


(3.75) 


and consequently the corresponding scalar potential Virf>, x) reads 


V(<f>, X ) = 


6 0 a(l + a)it s - <p) 1+a + \/ fc oa(l + a)(t s - (p)^ 1+ ") +b 0 a 


(1 + a)n 2 


+ —P5~{ts ~ X) 


U-l+a) 


2 

-6 0 (1 + a)(t s - (p) a - 2y/boa(t s 

-4>) 2 

4^ait s - x) 1+2 " } 


(1 + a)n 2 


(1 + a)n 2 


(3.76) 


Therefore, the two-field scalar-tensor theory that generates the cosmological evolution of 
(3.72), is given by expressions (3.74) and (3.76). 


3.5 The case of multiple scalar fields 

We close the current section by presenting, for completeness, the possible generalization of 
the above results in the case of multiple fields. Such a framework offers huge capabilities 
and possibilities of cosmological evolutions. We provide here the general framework of this 
extension following [106] (see also [128]). The scalar-tensor action in the case of n real 
scalar fields (pi , i = 1 ,..n, reads 


5 = 


V i=l 


ui(pi)d^cpid^(pi -Vicpi 


(3.77) 
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Assuming a flat FRW background, and assuming that the scalar fields are functions only 
of the cosmic time t, the corresponding Friedmann equations write as 

n 

V 'u((j)i)<jr = - 2 #, (3.78) 

i =1 K 

= -2 (3i7 2 ) - - . (3.79) 

1=1 


We choose the function f(t ) in order to satisfy the following equation: 

n 

y u;;(t) = -^/'(t), 

1=1 


and therefore we obtain 


Wj(0j) = 


2 d/(0i,<fe,-,<ftn) 

k 2 d (j>i 


V((/>!, (f)2, ■■■An) 


K 


3/(01, 4 > 2 , 0n) 2 + y 


(01 ? 02 , - - ■, 


1=1 


d^i 


(3.80) 


(3.81) 


with the function f(U) = ...,1) = /(f). The solution of the reconstruction method is 

then given by 

<t> = X = t, H = (3.82) 

By appropriately choosing some arbitrary functions g n ((pi) we can make the kinetic func¬ 
tions wj(0i) to be 

wi(0i) = —-2 {/'(0i) + 52 OM H-1- 5n(0l)} , 

HZ 

2 

<^2(02) = —2 ^2(02), 


^nAn) — 9 9n{,4 > n) • (3.83) 

KT 

Hence, one may choose to have ghost and non-ghost scalar helds and describe the cosmo¬ 
logical evolution in a desired way. We do not analyze this case in more details since it lies 
beyond the scope of the present work. 


4 F(R )~gravity analysis 

As a last section of this work, and for completeness, we will investigate the case of F(R ) 
gravity [101, 104, 105, 129-134], which is closely related to scalar-field theory, focusing 
again on the Type IV realization of a Hubble rate of the form (2.2). We will follow the 
reconstruction procedure of [135-137] (see also [69]). 

We start by recalling that in the case where a 1, the scalar-tensor potential which 
perfectly describes this cosmology is a hilltop-like potential of the form (2.13). We shall 


- 27 - 






assume that the parameters cq and bo appearing in (2.14) are chosen in order to satisfy 
c 0 < 1 and bo <C 1, and in addition a 3> 1. Furthermore, we remind that for an FRW 
metric, the Ricci scalar writes as 


R = 6(2 H 2 + H). 


(4.1) 


As the cosmic time t approaches the finite time t s , then t — t s —>• 0. We shall use this 
property later on in this section, in order to find approximate expressions for the F(R ) 
gravity. In addition, as t —» t s , the Ricci scalar approaches the limiting value R = 12Cg, as 
can be seen by calculating it’s analytic form 

R = 12cq + 246qCo (—t + t s ) a + 126q(— t + t s ) 2a — 6bo(—t + t s ) ^ a a . (4-2) 


We shall further investigate the physical consequences of our assumptions in the end of 
this section. 

The Jordan frame action of a pure F(R) gravity is 

S = J^J d A x^gF(R ), (4.3) 


with pure referring to the absence of matter fluids. Varying action (4.3) with respect to 
the metric tensor leads to the Friedmann equation 


- 18 


4 H(tfH(t) + H(t)H(t) F"(R) + 3 H 2 [t) + H(t) F'(R)- 


F(R) 


= 0 . 


(4.4) 


By making use of an auxiliary scalar field (which has no kinetic term), the F(R) gravity 
action of (4.3) can be rewritten as 


S = j d A xV=g[P{(t>)R + Qm ■ 


(4.5) 


By finding the exact analytic form of the functions P((j>) and Q(4>), which have an explicit 
dependence on the scalar field <f>, we may easily obtain the specific F(R) gravity. This can 
be done by varying the action of (4.5) with respect to the auxiliary scalar degree of freedom 
4>. Such a variation gives 

P\ ( j ) )R + Q'{ct>) = 0, (4.6) 

with the prime now denoting differentiation with respect to </>. Solving Eq. (4.6) with 
respect to 4> will provide us with the function 4>(R), and then the F(R) gravity can be easily 
found by substituting 4>(R) to the auxiliary F(R) action of (4.5). Hence, the reconstructed 
F(R) gravity will be 

F(<KR)) = P{<KR))R + Q(<KR)) ■ (4-7) 

Let us apply the above reconstruction procedure for the case where the Hubble rate is 


H{t) = c 0 + b 0 (-t + t s ) a ■ 


(4.8) 
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By varying action (4.5) with respect to the metric, we obtain the equations 


0 = - 6H 2 P(<p(t )) - Q(0(()) - 6 = 0 , 

0 = (4H + 6H 2 ) Pm)) + QWS) + + ^r 1 = 0 • ( 4 . 9 ) 

Eliminating Q((p(t)) from (4.9) results to 

2 d^«t))_ dpp) +4i ^ 

d t z dt 

In this way, given the Hubble rate, the explicit form of P(t) can be found and from this 
we can obtain Q[t) accordingly. We mention that, due to fact that the F(R) action (4.3) 
is mathematically equivalent to the auxiliary F(R ) action (4.5), the auxiliary scalar held 
is identical to the cosmic time t, that is cf> = t. This was proven in detail in the Appendix 
A of Ref. [135]. Substituting the Hubble rate (4.8) in (4.10), we obtain the following 
second-order differential equation: 

- 2 [ fi o + M-i + tsT] - 4&o(-* + t s )~ l+a aP(t) = 0 , (4.11) 

which by setting x = — t + t s can be rewritten as 

- 2 (co + b 0 x a ) - 2b 0 x~ 1+a aP(x) = 0 . (4.12) 

The general solution of the above equation is difficult to be extracted, thus we will follow 
approximations. In particular, since x is small the differential equation becomes 


0 d 2 P(x) 

dx 2 


2 ClJ ^=0. 

dx 


(4.13) 


with solution 

P{x) = -— J - + C 2 , (4.14) 

1 + Co 

with ci, C 2 being arbitrary constants. Then, the function Q(x) can easily be calculated 
using (4.9) as 


Q(x) 


-6 CqC 2 


12 cQCie C03: 
1 + Cq 


126 0 coC2X a 


186 oCoCie co 3 : x“ 
1 + Cq 


66 qC 2 x 2q 


CQXry.’i.a 


66qC] em~x 
1 + Cq 


(4.15) 


Substituting the hnal forms of the functions P(x) and Q(x) into (4.6), and by solving with 
respect to x, we can find the function x(R). For the case t —>• t s that we are interested in, 
i.e. for x —> 0, we can expand the involved P'{x) and Q'(x). Doing so we acquire 


p/ M „ c QCi , c p c i^ , CpCi^ 2 
1 1 + Co + 1 + c 0 + 2(1 + C 0 ) ’ 


(4.16) 
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and 


Q'(x) ~x 


2 a 


+ x v 


+ x u 


6 (6o c oCi) 

6(&o c o c i) x 

+ x a 

12 (6 0 CoCi) 

12 (6 0 cgci) x 

1 + Co 

1 + Co 

1 + Co 

1 + Co 

6 (frococi) 

6(frococi) 

X 

1 

12 (cg Cl ) 

12 (cqCi) x 


1 + Co 

1 + Co 


\ 

1 + Co 

1 + Co 


6(b 0 c 0 cia) 
(1 + c 0 )x 

6 (6 0 cgcia) 

1 + Cq 






+ x 


2a 


r 12 


(l+c 0 +C2 ) 

a 

12 (blcocia) 1 

l 

X 

1 + Cq 


12 


+ /<- 


boco 


c 1 

1+CO 


+ C 2 


a 


x 


12 (b 0 clcia) 
1 + Co 


Since a>l, keeping the most dominant terms in (4.17) we obtain 

Q'( x ) = (hi?* _ 124£i2\ 
v ; V 1 + CO 1 + CO ) 

Hence, substituting (4.16) and (4.18) into (4.6) we acquire 


x 


12 cq — cqR + \J 144c® — CqR 2 


c l R 


(4.17) 


(4.18) 


(4.19) 


Thus, the F(R ) gravity form that can reproduce the Type IV singularity, can be obtained 
straightforwardly upon substitution of (4.19) into (4.7), and reads as 


Ev m ^ 144 cqCi _ 144cqCi „ 12cpCiy/l44:CQ Cq R 2 _ 12cqCi y/l44c® CqR 2 

{> ~(l + c 0 )R (l + c 0 )R 2 (l + c 0 )i? (l + c 0 )R 

(4.20) 

Observing the above form, we deduce that by appropriately choosing the arbitrary constant 
c 2 we can make the F(R) function to be the usual Einstein-Hilbert gravity plus curvature 
corrections. Indeed, for c 2 = 1, the F(R) form becomes 


F( m ~ /? + 144c o Cl _ 144c o c i , 12c 0 ci \/144cq — CqT? 2 _ 12cgci Vl44cg - c 2 0 R 2 
1 j (i + co)^ (l + co)R (1 + co )R (l + c 0 )R 

(4.21) 

We can further simplify the final form of the F(R) gravity by recalling that (4.21) is valid 
only near the Type IV singularity, in which case the Ricci scalar approaches the limiting 
value R = R c = 12cq (see (4.2)). Therefore, for R —>• R c , the last two terms of (4.21) are 
nearly zero. Consequently, the final form of the F(R) gravity is approximately 


F(fl)^fl+ 144 f 1 cl(1 : co) i 

(1 + c 0 ) R 


(4.22) 


This type of F(R) gravity was described in [101]. Recall our initial assumptions, namely 
c 0 < 1 and 6 q T an d in addition a>l, Furthermore, as t —>• t s the scalar curvature 
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approaches the limiting value R —» R c . Therefore, we have an appealing physical picture if 
we additionally require that t s is a future late-time singularity. Indeed, the F(R) gravity 
is of the form F(R) ~ R+ ^ and the scalar curvature approaches a late-time point, which 
since Co 1 this limiting curvature value can be chosen to have a very small value. Hence, 
the F(R) form of (4.22) can successfully describe late-time acceleration near the future 
Type IV singularity, which was the goal of the present section. 

We mention that the F(R) gravity of the form (4.22) can consistently describe dark 
energy and late-time acceleration, as was demonstrated in [101], and with the presence of 
an R 2 term it can unify inflation with dark energy [138]. However, the model (4.22) is not 
completely in agreement with observations, having problems with the correct description of 
the gravitational interaction in the Solar System. Nevertheless, the model (4.22) turns out 
to be a very useful approximate form of consistent and realistic versions of F(R) gravity 
(see [133] for a review). However, we have to note that the de Sitter solution for this model 
is Rds = 12 cq - v /— 3 ci , thus for negative ci the model we found may not describe dark energy, 
as was shown in [139, 140], since F"(R) < 0. Additionally, as it was stressed in [141], in 
order for the scalaron mass to be less than the Planck mass, the F"{R) has to be bounded 
for a finite R, which is not satisfied by (4.20). 

5 Conclusions 

In this work we demonstrated that finite time singularities of Type IV, can consistently be 
incorporated in the Universe’s cosmological evolution, either appearing in the inflationary 
era, or in the late-time regime. We thoroughly discussed the most consistent framework 
in which these can be realized without instabilities. In particular, using only one scalar 
field then the linear-level instabilities can in principle occur at the cosmological time at 
which the phantom divide is crossed during the dynamical evolution, as we showed to 
happen in the case of a hilltop-like potential. On the other hand, when two fields are 
involved, one being canonical and one being ghost, it is possible to avoid such instabilities 
during the phantom-divide crossing. In addition, despite their increased level of complexity, 
the two-field scalar-tensor theories prove to be able to offer a plethora of possible viable 
cosmological scenarios, at which various types of cosmological singularities can be realized. 
For instance, it is possible to describe phantom inflation with the appearance of a Type IV 
singularity, and phantom late-time acceleration which ends in a Big Rip singularity. 

For completeness, we also presented the Type IV realization in the context of F(R) 
gravity, which is known to be connected with a scalar-field theory through conformal trans¬ 
formations. In particular, we showed that specifically reconstructed F(R) forms can gener¬ 
ate cosmological scenarios with a dynamical evolution corresponding to hilltop-like poten¬ 
tials. As we evinced, the resulting F(R) gravity is of the form R + and this can be quite 
physically appealing if the Type IV singularity occurs at late times, since the resulting 
model is known to produce late-time acceleration. 

In the absence of a complete understanding of the nature of cosmological singulari¬ 
ties, we believe that with the present study we demonstrated that not all singularities can 
be harmful, at least for the cases where observational consistency is the main goal of a 
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self-consistent physical theory. Indeed, even in observable quantities that contain higher 
derivatives of the scale factor, the effect of singularities, and particularly of Type IV singu¬ 
larities, might be negligible, as for example in the jerk parameter. Definitely, a lot of effort 
is needed to theoretically understand the singularities, but probably these Type IV singu¬ 
larities might be a classical bridge between classical cosmology and the ultimate quantum 
gravitational cosmology, which offers a complete description or at least a remedy of these 
classically unwanted spacetime points. Some steps towards these directions of research are 
offered by loop quantum gravity, in which no singularities occur [142-146]. However, we 
mention that it is possible even at the classical level to remedy singularities [23] or even 
avoid them [60], when quantum corrections are taken into account, as in [6, 23], and thus 
we hope to address this issue in the future in more details. 

It would be interesting to notice that phantom-non-phantom transitions as well as non¬ 
singular evolution appear also in the case of bouncing cosmology. Hence, one could follow 
the procedures and the techniques of the present work, in order to investigate the effects 
of Type IV singularities, and probably of Type III ones, in the context of various ekpyrotic 
[147-149] or bounce Universe scenarios [29-50]. This thorough analysis lies beyond the 
scope of the present work and will be addressed elsewhere. 

Finally, we have to mention that in this work we did not address in general the issue 
of possible analyticity loss beyond p = 0. As it is also commented by Barrow and Graham 
in [74], for large-field inflation potentials of the form ~ p n , there are fractional powers of 
n for which the potential is always real, and extensibility through p = 0 can be achieved 
leading to a Big Crunch (for instance when the potential becomes negative). The forms 
of the potentials we used are exactly of this class (large-field inflation ones), or hilltop-like 
potentials. Nevertheless, in general, the analyticity issue around p = 0 is very interesting, 
both from mathematical as well as physical point of view. This issue however deserves a 
detailed investigation of its own and is left for a future project. 
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